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ABSTBACT 


Thlc report docuBonc* and sunurlzes the aecoopllshaieacs 
over cho past year in two areee: (1) developnenc of Landsat claasl- 

fiction accuracy aasessoant techniques, and (2) development of a 
computerized system for assessing wildlife habitat from land cover 
maps. This report Includes a literature review on accuracy 
assessment techniques, a complete explanation for the techniques 
developed under both projects, including example analyses and 
listings of the computer programs. 

A summary of the presentations and discussions at the 
tfatlonal Working Conference on Landsat Classification Accuracy is 
included. Also, two symposium papers which have been published 
on the results of this project are included as appendices. 



1«0 Introduction 


Knny stud, m hnrn bnsn conducted to detotnlns ths usefulness 

of LANDSAT decs for sMpplng lend cover. However, very little reseereh hes 

been done to determine Che degree of success (i.e., eceurecy) in doing, this. 

A recent literature review by Meed (1977) indicated chat: 

...more work is needed to develop reliable techniques for 
esclmaclog classificscloo accuracies. A means of compsring 
Che accuracies (i.e., to compare classification matrices) 
obtained In different areas on different daces, or estimated 
by different techniques is needed. Such techniques should permit the 
investigator to test hypotheses that at specified level of 
confidence the accuracies from several areas, daces, etc. 
are not different. 


(p. 59) 

Mead (1977) continues by suggesting '^Future studies might consider 
iterative proportional fitting of the classification matrices as a 
means of doing this.” (Bishop ^ al . 1975). 

The apparent absence of quancicacive methods for comparing classi- 
fication accuracy is certainly a scumbling block chat must be overcame. 

The effects of imaging dace, spectral band combination, classification 
algorithm, training set selection procedure, and the image analyst on final 
classification accuracy must be studied. Therefore, the following study 
was proposed with these objectives: 

1.1 Objectives 

1. To develop a computer system that implements an iterative 

proportional fitting technique to "normalize'* the coefficients 
within classification error matrices. 
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2. To dtvolop hlorarehol aod«l« for CMtlag tbo •IgalfleoBco 
of Mvoral foetoro (o.g., iaago d«tt, elaaoif let cion 
algorlehm* tho analyse « tec.) on eht rtaulting claaalfl* 
eaeloa accuracy. 

3. To ttat eht about etchnlquta and dtetealnt ebtir uatfulntaa 
wleh aceual daea for claaaiflcaeion accuracy. 

1.2 Juaelflcaelon 

Sltatacch will undoubetdly concinut coward dtvtlopntnc of a ays eta 
for classlflcaelon of land eovtr from dlglcally rteordtd Landaae ioagtry. 
Such rtstarch tf fores will in part bt mtasurtd by Improvtmtnes in eht 
classlflcaelon accuraclts achitutd. Thtrafcrt scltntlsts will nttd ways 
of asstssing eht accuracy. Also eht accuracy of eht final maps produced 
muse bt utrlfitd btfort ehty art dlscribuctd eo users. Once standards 
are tscabllshtd« rigorous seacistlcal procedures will be needed eo aalnealn 
the quality of ehe maps. Therefore, ie can be seen ebae accuracy assess* 
menc techniques will be needed in both the research and operational 
environments. 

1 . 3 State of the Art of Landsat Classification Accuracy Assessment 
Landsat, like any ocher remote sensing system, is only as good 

as our ability to evaluate it. The need for techniques eo assess the 
accuracy of the Landsat sensor systems cannot be understated. As 
Freese (1960) states, "testing the accuracy of some measurement against 
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«n «ec«pt«d ttaadArd rtquirM * «c«t«m«ac of eh« oeeuroey roqulrod, * 
BMOuro of ch« oeeuroey octaiaod, and an objaeelva aaeboo of daciding 
uhachar eha accuracy aecalnad la equal eo cha accuracy raqulrad**. If 

0 

tiiara ara no aachoda for aaaaurlng cha accuracy aecalnad wlch a carcaln 
aanaor syacan, chan chara uiU ba no uay co aaka cooparlaona bacuaan 
ayaeaaa co dacamina which la bacear. 

If Landsac la evar co bacona an oparaclonal syseem, chan 
avaluaclon and accuracy asaaasmanc cachnlquaa oniac ba davalopad co 
show whara such sansoc syscans glva mora adaquaca rasulcs chan con- 
vanelonal machods. Thasa aaaasaaanc cachnlquaa auac chan ba appllad 
eo spaclfic appllcaelona. For axaspla, *%ha usafulnaas of saealllea 
Inagary for forascry dapanda on cha axcane eo which foraac daca can 
ba cacordad by a raaoca aanalng syacam froa saealllea aleleudas, pzo- 
cassad by an iaaga Inearpraeaelon syseaa* and oaad in forase oapplng 
and Invaneorlas" (Kalanaky and Schark, 1975). 

1.31 Accuracy Asaassnanc Tachnlcuas 

Thara hava baan vary faw acudlaa dona on accuracy of Landsac 
classlficaclon. Hose of cha aarly asaasaoancs wara dona as an "afear 
choughc" wlchouc much conaidaraclon glvan eo cha scaelsclcal oaehods 
usad. Thasa scudlas, such as eha ona dona by Kalanaky and Schark (1975), 
usually dealt only with training sat accuracy. Tha use of training secs 
as wall as ocher possible areas to ba assasaad will ba discussed lacar. 
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A review of the current eeeeeeBent techniques ere neceesery before 
eny of the epplicntions of these techniques esn be understood. 

The oosc conoon wsy to describe the sceurecy of e Lsndset iasge 
is in the form of en error aecrix (e.g., Todd ec , 1980; Heed and' 
Meyer, 1977; Hoffer, 1975). An error nscrix is a square array of 
nuabers see out in rows and coluans which express the nuaber of pixels 
assigned as a particular land cover type relative to the actual land 
cover as verified in the field or froa photos. The coluaos usually 
represent the grouiul truth and the rows indicate the computer assigned 
land cover category. This form of expressing accuracy 
as an error matrix allows for an effective way to evaluate 
both errors of inclusion (commission errors) and errors of exclusion 
(omission errors) present la the classification. Also, the error matrix 
allows the analyse to dataraine the performance for individual categories 
as well as for the overall classification (Boffer and Fleming, 1978). 

In the ideal situation, all the oon-aajor diagonal elements of the 
error matrix would be rero, indicating that no pixel had been aisclassi- 
fled (Lillesand and Kiefer, 1979). 

There are two basic types of accuracy assessments. They are site 
specific accuracy and non-site specific accuracy. All the methods 
described to assess accuracy can be applied to either type. Non-site 
specific accuracy is less useful chan sice specific accuracy. Meyer 
et al. (1975) used a non-sice specific accuracy assessment to evaluate 
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clM«ifie«tlon of Undsot ia«g«>nr la SottChaMCtirn Hontaaa. Xoeol aroa 
acraagat wara ealeulaead for aaeh laforaaclonal elaaa. Thara wara ao 
eaaca t-tda for poalclonal accuracy (alea spaciflc), Juac ralaciva eoeal 
acraagaa. Hayar found eha aaeiaaca of eha ralaciva proportion of aach 
covar cypa cotnparad favorably with tha ground truth (i.a. • actual acraa 
of aach land covar caeagory). Hoi^ar, ha alao ooclcad chat onlaaion 
and conaiaalon arrora wara vary obvloua and chat tha ovarall poalclonal 
accuracy of eha covar eypaa wlehln eha araaa aeudlad waa poor. 

Thla axaspla polnea out eha oajor dlaadvaneagaa of a non^alca 
specific accuracy aaaaaasane. If only eoeal acreage aatlmaeaa are needed , 
than chla method may apply. However » the natural raaourca manager la 
uaually Inearaaead In eha location aa wall as the acreage of a certain 
land covar caeagory. If ebla la the case, le la obvloua chat non-slea 
specific accuracy aasaasaane la noe adequate. 

Sica specific accuracy • on eha ocher hand* Is a aaaaura of how 
wall eha compucar (claaslfiu«eion algorleha) claaalflaa aach pixel with 
raspace to eha ground truth. It la a more meaningful rapraaancaclon of 
the accuracy of eha class If Icaclon. Tha analyst can sea which catagorlaa 
are easily Idenclflabla and which ar% being confused. Although Lyon (1979) 
used alee specific accuracy aaaaaatsanc , ha Includes no error macrlcaa In 
hla paper. Instead t ha gives juac one number as a measure of the accuracy. 
This Is a common problem throughout eha Ilearaeura. Without error matrices. 
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readsr has Uctla cha&ea of uodarstand ing how an accuracy figure 
was daeerminad. Tha raadar also losas Cha knowladga of which caca« 
gorias wara easily idancifiad and which wara difficult. 

Onca cha error aaerix has been generated, a very sisple procedure 
can be used to detervina cha overall accuracy. Since all the values on 
Che major diagonal raprasanc chose pixels chat have been correccly 
classified, if one adds up cha major diagonal and divides this number 
by Che cocal number of pixels classified, one will obtain che overall 
accuracy of thac error matrix. This is che most common use of che error 
matrix in accuracy assessmenc. 

In recent years, some new techniques have been developed to assess 
classiflcaclon accuracy. Among these new methods are analysis of vari- 
ance techniques, regression analysis techniques, and discrete multi- 
variate analysis techniques. Each of these methods has certain assumpt- 
ions chat must be met before the technique can be used for assessing 
classification accuracy. If these assumptions are not net, the technique 
loses its power. 

The data used in classification accuracy assessment is of che 
discrete type. Discrete data, as opposed to continuous data, may take 
on only a limited number of distinct values (Snedecor and Cochran, 1976). 
In analysis of variance, the data must be normally distributed in order 
to meet che assumptions of che technique. Since discrete data is not 
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oormally diseribue«d, le uouid t««m that AHOVA 1« aoe a good eachniqua 
for accuracy assaasma&t. Howavar. Rosaafiald (1980) haa propoaad cha 
ttsa of cha logit eranaforaaelon or cha aresina eransformaclon as daaeribad 
by Snadacor and Cochran (1976) to transform cha data into an approxlnataly 
normal dlacributlon. Soaanfiald atataa, "tha statistically incarpratad 
results of cha waighcad adjustmant agraa fairly well with what might ba 
cachnologically axpaccad, and ara charafora judged cachnically accept* 
able". After cha transformation is applied to the data, cha analysis 
of variance can ba run. From cha resulting AliOVA cable, multiple range 
eases ara applied eo population means found to ba significantly different 
(Rosenflald, 1980). Analysis of variance is a powerful staclstical cool. 
However, ocher cechnlques chat do not reqtilre so much data manipulaclon 
should also be tested. Rosenflald (1978) agrees, ''this does not mean 
that they (ANOVA) are Che best; however, Che cools available should be 
used until something becter comes along". 

Regression analysis is another way of vlsiially representing 
accuracy. In this case the ground truth (i.e., actual laxid cover) is 
the independent variable, X, and the computer classification is the 
dependent variable, 7. If the computer is completely correct in its 
classification, then all the points will lie on a forty five degree line. 
More likely, the points will be spread out from this line. The value of 
Che correlation coefficient can Chen be used Co get an idea of che 
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r«l«tlve agr««a«ne b«cv««n th« ground truth and tha conputar claaai* 
fication. Ragraaslon analysis has not baan wldaly usad in tha lltara- 
tura and tharafora no mora will ba said about it. 

1.32 Sampling Taehnli^tiM 

The need to use nora than just training areas for accuracy 
assessment has already been discussed. Howenrer, one could not afford 
nor desire to assess the entire seme. Instead, a representative sample 
should be chosen and assessed as the accuracy for the entire scene. 
Sampling allows not only the calculation of a number that represents 
the accuracy of the classification, but also allows for a confidence 
interval to be placed around that number. 

Ginevan (1979) states three criteria that should be satisfied in 
any sampling scheme. -These criteria are: (1) the sampling scheme should 

have a low probability of accepting a map of low accuracy, (2) the sampl- 
ing scheme should have a high probability of accepting a map of high 
accuracy, and (3) the sampling scheme should have a minimum number, N, 
of ground truth samples. }lany researchers (Hay, 1979; Ginevan, 1979; and 
Genderen and Lock, 1977) agree that stratified random sampling is the 
best sampling scheme to use. Rhode (1978) proposes ocher schemes 
including cluster-stratified sampling and two phase sampling. No matter 
which sampling scheme is used, it should be chosen so to obtain the 
maximum information with the minimum amount of work. This Involves 
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eonsid*rlng oany varlablM such at tarraln, iaaga Identifiabla 
loaeloaa* and variability of land covar caeagorlaa. 

It should also ba noted that errors arise in classification 
from ocher sources besides the sanpling scheme chosen. Problems arise 
in radiometric correction and geometric rectification. Also, the time 
interval betwem when the Imagery is attained and when the field check- 
lag is done may cause differences in land cover category. It must also 
be reallaed chat just because the classification of a category seems 
perfect, this does not always mean chat the method is error free. The 
result may occur purely by chance because of the sampling design. "This 
fact is seldom appreciated by many image interpreters when checking the 
accuracy results of their remote sensing land use survey (Genderen and 
Lock, 1977). 

Finally, ho matter which sampling scheme is chosen, a sample 
size must be determined. This situation is described by Glnevan (1979), 
'The sampling problem as defined here is the determination of the optimal 
number, N, of ground truth samples and an allowable number, X, of mis- 
classifications of these samples." Once these have been determined, the 
results of image interpretation are checked against the N ground truth 
samples and the map is accepted as accurate if X or fewer of the ground 
truth samples are mlsclassified. The optimum number of samples, N, to be 
taken has met with widespread disagreement throughout the literature 
(Todd et al. . 1980; Hay, 1979; Genderen at al. . 1978; Genderen and Lock, 
1977; and Hord and Brooner, 1976). Each researcher seems to have his 
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o«n Ideas about sanpla also deeetainacloii and it la obvious chat a 
great deal more research Is needed in this area. 

1.33 national Data Base for Error Hatrices 

Letters were sent out to potential sources of error oatrices 
asking that any fflatrices they had be sent to us for inclusion in a 
national Data Base for Error Matrices. An infomatlon questionnaire was 
sent along with each request for data. This questionnaire contained 
questions about the location of the area analysed, the analyst, the 
algorltha, and the data the data were taken. 

All error natriees chat we have received have been compiled 
along with their corresponding pertinent information and placed on a 
computer cape. This data are available for distribution to ocher users 
upon request. A listing of the sources of error matrices can be 
found in Appendis I. 


Two tMthod* w«r« us«d co coop«r« cwo*diaenslonAl maericos 
raptMaaclng ground elMtlflencloa vorsus machine classlficacion 
from dlffaranc methods. In the first method of comparison the 
cell entries in each matrix are successively balanced until the 
sum of each of the matrix margins is one. The entries in the 
matrix then represent a normalized percentage of the total 
observations occurring in each matrix cell. Within an individual 
matrix these percentages can be used to examine omission and 
commission errors. Classification errors between two or more 
machine classification methods can be evaluated by comparing the 
percentages in corresponding cells in each matrix. Matrices with 
differing numbers of observations can be compared since the entries 
in each matrix are transformed to percentages. 

The second method of comparison was a measure of agreement 
for two-dimensional square matrices presented by Bishop ^ (1975) 

A 

This measure, R, is calculated as the difference between the actual 
agreement and chance agreement between two classification methods. 

In this application the two methods are ground classification and 
machine classification. The measure is calculated as 
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. NEX - £ X X 

X - 1«1 1«1 

wh«r« r Is the aud>sr of rows in ehs matrix, X is chs number of 

ii 

observations in row 1 and column i, X^^ and X^^ are the marginal 
totals of row 1 and column 1, respectively, and M is the total 
number of observations. An approximate large sample variance, 
based on the asymptotic normality of X, is available, and can be 

A 

used to derive a confidence interval for X from a single matrix 

A 

and to perform tests for equality of X between two matrices. 

The two methods described above can be used together. Method 
two, X, will indicate whether two matrices exhibit the same degree 
of classification success Cor error). If a difference exists, 
method one can be used to determine in which particular category 
or categories the difference lies. 

2. 2 Categorical Data Analysis 

The influence of factors such as season of imaging, film type, 
and interpreter bias on classification accuracy was examined using 
categorical data analysis (Bishop ^ al. , 1973). Using this analysis 
technique the dependence of classification accuracy on a single 
factor or combination of factors cui be assessed. 
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C«c«gorlcal data analysis raqulras only that aaeh factor 
balng sxamlned for influanca on accuracy can ba assigned to an 
unambiguous category within each factor. These categories may 
be normative, ordinal, or interval. The result of data 
collection la a multidlmenaional matrix with each factor. Including 
ground and machine classification, serving as a dimension of the 
matrix. 

This method of analysis avoids the more restrictive 
assumptions inherent In alternative analysis methods such as 
multivariate regression or analysis of variance. !fo normality 
assumption Is necessary, no factors need be considered as con- 
tinuous, and interprecaeion of many duasiy variables Is avoided. 
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3.0 SaatPl* Data Aiulys— 

3.1 MAaGFIT AoMlyU 

* 

A» previously discussed, the FORXRAM eooputer progrsa MAR6FIT 
Csee Appendix tl) iaiplements e notnslisecion procedure which stendardires 
each error oatrix for purposes of comparison. The accuracy of the 
classificaclon can chan be represented as a normalized overall perform- 
ance. This value Is calcttlated the same way as in overall performance 
(i.e. . summing the major diagonal and dividing by the total) except 
that the matrix is normalized first. 

Smith and Itkowsky (1978) compiled five error matrices for a 
study in north central Colorado. Two of the matrices were for 
training sets; Original was compiled using a supervised classification 
idille Josesigs was compiled using a modified unsupervised classification. 
The ocher three matrices (Scrambll, Scrambl2, Scrambl3) were attempts 
to reclassify incorrect pixels using a computer program called SCRAHBL. 
Table 1 shows the Josesigs error matrix before normalization and 
Table 2 shows the matrix after normalization. 


Computer Classification Computer Classification 



Tablft 2. JosMlgs «rror saerlx after aormallzaelon. 


itaference Data 


1 

J 

Decld. 

Coolf . 

Grass 

Meadow 

Shrub 

Water 

Sa^e 

1 

Decld. 

.5363 

.0996 

.1001 

.0168 

.0645 

.0155 

.1674 

Conlf. 

.1382 

.8044 

.0158 

.0027 

.0102 

.0025 

.0265 

Grass 

.0154 

.0200 

.5025 

.0842 

.1943 

.0156 

.1681 

Meadow 

.0879 

.0035 

.0174 

.7139 

.1457 

.0027 

.0291 

Shrub 

.1804 

.0180 

.0906 

.1366 

.4089 

.0141 

.1515 

Water 

.0035 

.0046 

.0230 

.0039 

.0148 

.9108 

.0385 

Sage 

.0384 

.0499 

.2505 

.0420 

.1615 

1 

.0389 1 

.4190 , 


Normalized 

Overall 

Perfonoaace 


4.2958 


.61 


r«4.2958 
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lAbl* 3 shows ths rssules of ovsrsll psrforaanes and 
noriBsllssd ovsrsll psrforssacs for sll flvs srror nscrlcss. Noes 
chae chs rslstlvs seeursclss srs sisdlsr for ebs evo psrforasaes 
vsluss sxespc for chs Josssigs astrlx. Csrsful study of Tsbls 1 
shows why ehis is so. Only thrss plxsls In ehs shrub estsgory wars 
corrsecly clssslflsd. This forcsd chs norosllssclon proesdure to 
loflses chs vsluss In chs shrub row snd coluan dscrssslng chs 
nomsllzsd psrforosncs accuracy. Also, no sags caesgory plxsls wars 
clssslflsd at all rssulclng In cha sasa cyps of noroallzaclon 
problsm. 

Tabls 3. Ovsrall and norsalizsd ovsrsll psrforaanes rssules 
for flvs clasalflcaclon srror aaerlcss. 

Noraallzsd 


Hacrlx 

Ovsrall Psrforaanes 

Ovsrall Psrfozaanca 

Original 

90.371 

86.03X 

Josssigs 

8S.96X 

61.36Z 

Scrambll 

8S.43X 

79.97Z 

Scraabll 

78.94Z 

70.49Z 

Seraabl3 

80.18X 

74.17Z 



17 


SiallAt tMulet wtr« aehltv«d for aaetleoo eo^llod by Roffor 
C197So). tero four orrer maerleot «ur« eoapilod ae two difforoac dac«s 
eomparlftg a elaaaifleaeion of aajor land eovar eypoa varaua foraac 
covar cyyaa. Tha raaulea of ttoraallzaelon ahown la Tabla k agraa with 
tha ovarall parfomaaea valuaa ealeulacad by Boffar. 

Tabla 4. Ovarall aod oormaliaad ovarall parforaaoca raaulea 
for four covar eypa arror aacrlcaa. 

Itoroaliaad 


Matrix Ovarall Parforaaoca Ovarall Parforaaaca 


Major Land 
Covar Typaa 
6-5-73 

85.96Z 

89.51Z 

Major Land 
Covar Typaa 
8-8-73 

69.3SZ 

72.53Z 

Foraac 

Covar 

Typaa 

6-5-73 

71.79X 

76.875 

Foraac 

Covar 

Typaa 

8-8-73 

48.83Z 

57.885 


3.2 KAP?A Aoalyala 

Tha FORTBAM coapucar program KAPFA (aaa Appandlx III) calculacaa 
a K aeaelatlc for a glvaa arror matrix which allowa ooa to coapara 
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•rtor maerlcM to ••• If ehoy «r« •ignlfleancly difforonc. This typo 
of eofflporloon hM ooay uoot. I& on oxosplo oitod obovo, Hof for (197 5o) eooi- 

A 

pllod two eloaolfieotions ot evo difforone docoo. Tho R ocociotle and 
cotTo«?<mdlni coofldoaeo iatorvol (i.o.* uppor ond lowor bounds) oro 
prooontod for ooeb orror nocrix in TOblo 5. 

Toblo S. R seociotle vlch uppor ond lowor llaico oc 95S 
eonfldonco Incorvol for four covor typo orror 
moerleoo. 


Matrix Lowor Lialt R Dppor UjdIc 


Major Land 
Covor Typoa 
6-5-73 

.69396 

.69458 

.69521 

Major Land 
Covor Typos 
8-8-73 

.62880 

.62929 

.62978 

Forost 

Covor 

Typos 

6-5-73 

.38961 

.39055 

. 39150 

Forost 

Covor 

Typos 

8-5-73 

.33004 

.33074 

.33144 
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At ctn bt ttta from Ttblt 5* non* of eht confidtaet iaetrvtlt 
ovtrltp; chtrtfort, til ebttt mterlctt trt tlgnlfletiitly dlfftrtnc. 
Thlt mttat chte cht iatgery etkta te ewo dlfftrtac dtett it tignlfl- * 
cttely dlfftrtac which iaplitt chte oat dtet mutt chta bt bttttr ehta 
eht oehtr. A quick look te eht dtet iadletett ehtc 6~5~73 wtt cht 
•ignifictnely btcctr dtet. 

Anochtr txampLt of chit etchniqut it provided by Hofftr (1973b). 
In chit txtaplt, four mterlctt wtrt gtntrtctd from four dlfftrtnc 
cltttifictclon tlgorlehmt. Tht rttulcs prtttnctd in Ttblt 6 show 
chte til cht mterlctt trt tlgnlflctncly dlfftrtnc. 

Ttblt 6. K scteltcic wlch upper tnd lower limlct tc 932 


confidtaet inctrvtl for four eltttiflctelon 
tlgorlchms. 

* 

Mtcrlx Lower Lisd,e K Upper Limic 


Kontuptrvlstd 
(10 cl.) 

.60271 

.60479 

.60686 

Kontuptrvlttd 
(20 cl.) 

.38348 

.38373 

.38799 

Modified 

Supervised 

.47326 

.47381 

.47837 

.Modified 

Clusetr 

.71631 

( 

1 

.71846 

1 

. 72001 

1 1 


* 

A fintl txtaplt of cht K scteiselc it found in Appendix V. 
This txtmplt dttlt wleb compering phoco inctrprtctrt co set if chty 
trt tignifictnely different. 
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3.3 COWIABLE Analysis 

Thft APL computer program CONTABLE (see Appendix IV) allows 
one to analyze attltl>«ay contingency tables.* In the example here 
a 5-way table is analyzed. This table (Cameggie, 1972) deals with 5 
factors or effects listed in Table 7. The data consists of 18 
5x5 error matrices with various films, dates, and interpreters. 

Table 7. List of factors and effects for 5-way contingency 
cable. 


FACTOR 

EFFECT 


1 

Date 

(6/10, 7/25, 10/25) 

2 

Film 

(Color, CIR) 

3 

Interpreter 

(#1, #2, #3) 

4 

Row 

(1, 2, 3, 4, 5) 

5 

Column 

(1, 2, 3, 4, 5) 


The hypotheses to be tested in this example are listed in 
Table 8 while the results and conclusions are listed in Table 9. 


^Without the use of this program and its Iterative Proportional 
Fitting Procedure, analysis of tables larger than 3 dimensions 
would be Impossible. 
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Tablft 8. List of hypothosis for CONTABLE example* 



0 

0 

0 

^5 

0 

0 

0 


No film effecc 
No locerprecer effecc 
NO dace effecc 
0 No row-column effecc 

No dace-film inceraccion 
No dace-incecprecer inceraccion 
No filffi-incerprecer Inceraccion 


Table 9. 

Lise of 
example. 

results and conclusions 

for CONTABLE 


HYPOTHESIS 

CHI SQUARE VALUE 

CONCLUSION 


“o= “2 * 

0 

623.487 

reject Hq 


“o' “3 • 

0 

613.142 

reject Hq 


“o' “1 ■ 

0 

591.543 

reject Hq 


“o' “4"“j 

- 0 

134.485 

fall to reject 

«0 

“o' “12 • 

0 

145.961 

fall to reject 

«0 

“o' “13 * 

0 

162.393 

fail to reject 

«0 

“o' “23 ’ 

0 

144.707 

fall to reject 

«0 
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Table 9 shows chat alchough no single factor significantly 
affects the classification, the combination of two or more factors 
does. This means chat none of die three factors (film, date, 
interpreter) is more Important than the others. Instead all three 
factors Interact together to give the best classification. From 
Che analysis so far there is no significant one factor on which 
most of Che accuracy depends. 
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4.0 Accuracy ConfTenca 

A National Working Confaranca on Landsat Classification Accuracy 
Assassaant Procaduras was bald in Sioux Falls » South Dakota. A 
sussnary of this confaranca is givan in Appandlx VI as a draft manu- 
script which will ba ravlsad and submittad for publication in a 
Journal. 

5.0 Wildlifa Habitat Assassaant Mathods 

A sacondary task in this year's plan of work was to davalop 
digital spatial analysis techniques for assessing wildlife habitat. 
Appandlx VII includes a FORTRAN computer program for doing this, and 
Che techniques are described in a icanuscript which is Appendix VIII. 

6.0 Effects of Classification Accuracy on Interspersion Maps 
Artificial land cover type maps were made in order co test the 

effects of classification accuracy on computer generated interspersion 
maps. Three cellular maps were made, each containing 10 rows and 10 
columns with each cell assigned to one of S classes. The first map was 
used as a reference base map for comparison with the other two maps. 

The second map had 90Z of its cells classified similar co the first 
Ci.e., 90Z accurate), and Che third map was 70S similar co the reference 
base map. Five cover types, designated 1 through S, were used on 
each map. Simllarlcy matrices were generated between the accurate 
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(90S siallar) "Map II" and tha baaa map (Tabla 10), and bacvaan tha 
laaa accuraea, "Map III" C70Z aioilar) and eha baaa map (Tabla 11). 

KAPPA was uaad eo compara tha cuo rasulcing similarity matricaa. 

Tha Intarsperslon Indax daacribad by Maad ^ al. in Appandix VIII 
was uaad to craata Intarsparsion maps from each of tha three fictional 
cover type maps. The maps delineate areas of high Cdesignatad 3), 
medium (2), and low CD Intersperslon. Similarity matrices were 
created by comparing each of the Intersperslon maps (from the cover 
type maps II and III) with the Intersperslon map made from the base 
map (Tables 12 and 13). 

The implementation of tha KAPPA program (see Section 3.2) was 
than used to test for a significant difference between the Intersperslon 
maps. The resulting KBAT values indicate that cover type maps II and 
III were significantly different. A significant difference was also 
found between the two matrices for the Intersperslon maps. However, 
further work is needed to understand the effect of map accuracy on 
computer generated Intersperslon maps. Juxtaposition maps, and spatial 
diversity maps. Also, the effect of Increasing the number of cover 
types or the number of Intersperslon classes (high, medium and low) Is 


unknown 


Map III Classification Hap II Classification 
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Table 10. Slfflilaricy matrix for five fictional cover 
typaa on the baae map and on Hap ZZ. 

Baaa Map Claaaification 



1 

2 

3 

4 

5 

1 

19 



1 


2 


18 

1 



3 


1 

16 



4 

1 


2 

27 

1 

5 




3 

10 


Overall Accuracy "Jqq ■ 90* 


Table 11. Similarity matrix for five fictional cover type 
maps on the base map and on map ZZZ. 


Base Map Classification 



1 

2 

3 

4 

5 

1 

14 ' 

1 


2 


2 

2 

16 

4 

3 


3 

1 

2 

13 

2 


4 

3 


2 

17 

2 

5 




6 

10 


Overall Accuracy 


100 


702 























2 < 


Table 12. Slallarlcy matrix for three cacegorlea of 
Interaparsioi^ high (3) . medium (2) , and 
low (1) produced from the base map and 
map II. 


Base Hap Interspersloa 


EO 

i ^ 

« 
u 
0) 

^ 2 


Table 13. Similarity matrix for three categories of 
interspersion, high (3), medium (2). and 
low (1) produced from the base map and 
map III. 


Base Hap Interspersion 


U 

o. 

S 1 


^ 2 


i 3 
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7.0 Suamary and Futuf Watk 

Tha Iltaratura rtviav and prallalnary Invaaclgaciona show 
chat: CD eha ataclaclcal cachnlquas inieially proposad aca sound 

and aca usaful foe analysla of Landsac claasiflcaeion accuracy data, 
(2) subscanelal amounts of data from accuracy aasassmanca axisc 
but few sacs ara comparabla prohibiting hypochesas from baing 
castadf C3) prallmlnary rasulcs show chat tha machod used in sampling 
a classification can significantly affect the estimated accuracy. 

An ''automatic" computerised system needs to be developed for com* 
piling error matrices for any classification given the necessary 
ground truth and a specified sampling strategy. Experiments need 
CO be designed in the future so that fundamental questions can be 
answered about factors which affect classification accuracy. 

The wildlife habitat assessment system has greatest potential 
whan animals with requiraments related to the spatial characteristics 
of Che landscape are considered. Juxtaposition can be of great 
Importance or of very little importance depending upon the specific 
geographic area and the wildlife species of interest. When this data 
on Che spatial characteristics of the landscape are coupled with basic 
land cover information and ancillary data (e.g. , elevation, slope, 
soil type, political or ownership boundaries), an over*all system 
for habitat assessment may be realized. Such a system could be 
implemented on a computer and merged with data on other resource 
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accribuces («.g.. cinbar producing capablllcy). Fiirtbar work should 
Includa pilot tasting tha systam and an avaluation by fiald laval 


rasourca managars. 
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Appendix I. 

List of Sources of Error Hacrices 

1. Mead, Roy A., Landsat Digital Data Application to Forest 
Vegetation and Land-Use Classification in Minnesota. Ph.O. 
Disserts cion. University of Minnesota, 1977. 

4 aatrices (training set, test set, 2 evaluation areas) 
Minnesota. 

2. Fleming, Hiense) . Computer Aided Analysis Techniques for an 
Operational System to Map Forest Lands Utilizing Landsat MSS 
Data, LARS Technical Report 112277. 

2 matrices Colorado. 

3. Smith, James and Frank Itkowsky, Sensitivity of Variable 
Probability Sampling Estimates to Initial Landsat Classifi- 
cation, Final Report R.M.F. & R.E.S. USFS Coop-Agree. 16-741-CA, 
September 1978, CSU, Fort Collins, Colorado. 

5 matrices ( training sec, test set, 3 evaluation areas) 

Colorado. 

Madding, Robert and Harland Hogan, Detection and Mapping of 
Spruce Budworm Defoliation in Morchem Wisconsin Using Digital 
Analysis of Landsat Data. Proceedings of ASP Convention. 

Feb. 26 - Mar. 4, 1978. pp 285-300. 

2 matrices (normal and collapsed) 


Wisconsin. 
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5. Voss, A. W. , J. E. Bsksr, 6. E. Haussr, and 0. W. Nswton, The 
Uss of Lsodsac Dsrlvsd Land Covar Data in a Flood Peak Correla- 
tion Study, Proceedings ASP, Feb. 26-Mar. 4, 1978, pp. 135-146. 

2 matrices (normal and collapsed) 

North Carolina - Tennessee. 

6. Hoffer, Soger, Natural Resource Mapping in Mountainous Terrain 
by Computer Analysis of ERTS-1 Satellite Data, LARS Research 
Bulletin 919. Info. .Note 061575. 

10 matrices (different classification systems) 

Colorado. 

7. Hoffer, Roger, Computer-Aided Analysis of Skylab MSS Data in 
Mountainous Terrain for Land Use, Forestry, Water Resources, 
and Geologic Applications, LARS Info. Note 121275, 1975. 

4 matrices (varying spectral bands) 

Colorado . 

8. Hoffer, Roger, Mapping Vegetative Cover by Computer Aided 
Analysis of Satellite Data, LARS Technical Report 011178. 

2 matrices (test sites) 

Colorado. 

9. Hoffer, Roger, Variables in Automatic Classification over 
Extended Remote Sensing Test Sices, LARS Information Note 061571. 
1 matrix (test site) 


Indiana - Illinois. 
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10. Hof far, Rogar, Basie Poraat Covar Mapping Using Digicizad 
RaiBota Sansor Oats and ADP Tachniquas, LARS Information Hoca 
030573. 

13 matricas (casts at diffaranc spactral bands) 

11. Hallar, R. C. , R. C. Aldrich, R. S. Driscoll, R. E. Francis, 
and F. P. Wabar, Evaluation of ERTS>1 Data for Invantory of 
Forast and Rangaland and Dataction of Forast S Crass. PSff & RM 
For & Ranga Exp. Sea. Aug. 9, 1974. 

12 matricas. 

12. Ernst, Carols Lisatta, Digital Procassing of Remocely Sansad 
Data for Mapping Wacland Coontunlcies, Ph.D. Oissercation, 
Purdua University, Dec. 1979. 

6 matricas (classification) 

Indiana. 

13. Nelson, R. and R. Hoffar, Computer Aided Processing of Landsac 
MSS Data for Classification of Forest Lands, LARS Technical 
Report 102679, 1979. 

12 matrices 
Colorado. 

14. Carneggie, 0. M. , Large Scale 70 mm Aerial Photographs for 
Evaluating Ecological Coodieions, Vegetacional Changes, and 
Range Sice Potential. Ph.D. Dissertation, University of 
'.xlifomia, Berkeley. 

18 matrices (photo interpretation) 
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15. Lau«r. Donald, Clair a Hay and Andrew Banson, Quancitaclva 
Evaluation of Mulclband Photographic Techniques, Final Report 
for Earth Observation Division Manned Spacecraft Center, 

HASA Contract HAS 9-9S77, 1970. 

79 matrices (photo interpretation) 

16. Bryant, Emily and Gibb Dodge 
1 matrix. Maine. 

17. Roberts, Edwin 

1 matrix. Colorado (test set for bounty). 

18. Roller, Horman and Larry Vlsser, Acciiracy of Landsat Forest Cover 
Type Mapping in the Lake States Region of the U.S., Fourth 
International Symposium on Remote Sensing of Environment, 

April 23-30, 1980. 

1 matrix (Forest Cover Types) . 

Michigan. 

19. Newcomer, Jeffrey 

3 matrices. Pennsylvania. 

20. Harrington, John A. and Charles W. Dunn, Jr. 

3 matrices (forest - other) Oklahoma. 
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Appendix II. Listing of FORTRAN Computer Program MAR6FIT 


1 

3 

4 

s 


7 


9 
17 
I I 
!<> 
I i 

1 tt 
IS 

i 

1 7 
I ■» 

1 9 
■2'J 

as 

ay 
- > 


//‘naTFIV ,paGES «30 

c •«•••• ****•*••*••***•*••«•••••**••*•• ••***•••* 

c * * 

c • MAiiGFIT Ke^vklTTEN ANO OOCUME^^eo * 

C * wodSELL »i. CONGALTQN • 

C • OEPf. OF FORESTRY, VPlASl • 

C • JULY l«# 7 R * 

C * • 

c ***••••****•*♦*•*•••••*****•#**•**••••*•*****♦ 

c 

r 

c • 

c * This prusram .ias oesiuNEi) t(» change a 'aat^i* cf •-'axi-^um orAENSiu'^s 
C • uF S '»*50 I.'JTC A . 4 ATRIX .WTh nq^OE T?RmI NtO A nO Cni.'JN;N •'*a«GIn\l 3 

C * 

C * NuhmaT s The NUMnt.) OF haTBICFS 7 " h£ ChANGFC 

C • TA(iilI,J) s the value in Rij.-. I anC COL'j-"'* J OF I ^E GIVEN ■••Ar,<ix 

C • WHAPd) s THfc i»«ARGTMAL value fur ROia I 

C * C-tA.^(IJ a rnfe xAHGINAL VALUE FfR COLU <-»4 J 

C • '*A»ir a THE WAxIMsjM .NUWb£r? OF ITERATIUAiS 

C • MAXOEV a Th£ VAXl'Mir* ALLOaAhLE UEVIATIiIN 

c * IH a The NyWeER OF RO.»S Ii 4 THE haT.VIx 

C * IC s the Htjr-(*ER PF COLU»*'VS I-V Th£ , 4 ATqi» 

C * 

c ***•*••*•**••**•*•**•**•**••••**•***•*•**•»*•••***•*»**••******•*••• 

c 

c 

0 iHEMSrOf-i TAH (50,50) ,Ri^Aq(50) , C'**AR ( 5o ) , F I T { 5C ) 

REAL maxOEV 
C 

NCUUNTaO 

V£Ani 5 ,i 5 j 'IIJMV 4 T 

15 FQt,MAT(I£) 

C 

50 REA0(5» 10) H, IC 
10 FOkH^TCilS) 

ou aoo i8i,iR 

9 EA 0 ( 5 , 20 ) (TAHfI,J) ,jal ,K) 
a') FO<^"-A T(U(Fo, I ) ) 

a»» cO''inNu£ 

»kAu( 5 , 30 ) iRf’Akd), lal, £tr) 

30 FOkMAT( liCFS.O; ) 

qbAulS,aoj (C'^Ak C [) , lal , IC) 
ttO FO 3 >»AT( 12 CFN, 0 ) ) 

3 EA 0 ( 5 , 51 ) ■•'AX r T,‘tAXrEv 

51 F-J* vAr(IS,FI 0 . 3 ) 

9 KA 0 (o,S£) 

EOh^jTC'AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA.'AAAA'J 

L 

/•»! TE 1 ^, 99 “) 

99« F Jk .'A T { ' 1 ' ) 

^RITF. (b, 5 .i) 

.%-*IT£(S, 53 ) 

53 F(j 4 N-Ar (lx, '*•••***••*•*••••••♦•**'*••*****•••****«*'* *////) 

vRir£(S,Saj 

54 F0kM4 t ( 1 X , » TmE LHlOriAL ••■ATRT'f is: ) 

--7 I ft (N, 55 ) 

55 Fl... 1A f ( 1 X . ♦ _ _ _ *//) 

00 a 50 lai , 
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$0 


31 

56 

32 

250 


C 

33 


Stt 


35 


38 


37 

400 

36 

300 

19 

100 

aO 


ai 


'42 

5oO 

43 

150 

aa 


45 


46 


47 

700 

40 


49 


50 


5 1 

800 

52 

600 

53 


54 


55 


56 


57 


56 

1000 

59 


60 


61 


62 

900 

63 


64 


65 


66 

60 

67 


6A 

260 

69 

70 

70 


71 


T2 

60 

75 

1100 


c 

74 

1200 

75 

r 

76 


77 



Nirsu 

00 300 lsl,iw 

Fir(i)«o. 000001 

00 «00 jsific 


♦TA8(I,J) 


UU 500 J«i*ic 

rAfl(i,j)«rA8(i,j)*0MA»(i)/FiT{n 


t)0 600 JSUIC 
630,000001 
00 ^00 l*l»II< 


ou 000 i 8 i,iw 

IF( TA 6 ( I, J) ,Lr,l,oe-lO) TA6(I, J)s0.0 
TA3(I.J)3TA6(I,J}*CMA«(J)/6 


NlTaMT^l 

OsO.O 

DO 900 1st, in 
FIT(I)sO, 000001 
00 1000 J31,IC 


IF(0.6€.h) 60 TO 900 

OSM 


IF(0,L6*»«AXOEV) 60 TO 260 
IF(>»IT,l£.maxiT) 60 TO 100 
•<niT£(6f 60 ) **AXfT,U 


AIMUM oeviATlOM is: ',F10,3) 
60 TO 1200 


»I 5 »* I reftATICAS. */* THE CUWnEMT MAI 


AIATIOM of: SFIO.S/) 

00 1100 isi, m 

«»<ITe( 6 »« 0 ) (T46(I,J) , Jsl, IC) 


M 5 .* 


ITfcWATl(.^S /»TTH A MA<I-«UM U£\ 


IF (•''tCUONT.UT.MUwWA T) 60 To 60 


STOM 

END 


//DATA 


^!?ro. 




'?fc 




'^iTy 
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Appendix III. Listing of FORTRAN Computer Program KAPPA 


l 

a 

1 

a 

5 

ft 

7 

a 


9 

10 
1 1 


la 

13 

u 


IS 

Ui 

17 

la 

19 

ao 


n 


ai 

a« 

as 

ae 

37 

aa 

39 

30 

31 


//IfiaTFlV 

C 

C 


C 

C 


.P46ES33S 


* 


KAPPA WAS PEMPtrTEN ANO OQCUMENTEO 87 
RUSSELL 6, CONGALTQIM 
OEPT, OF FORESTRY, VPIASl 
JULY 1979 


C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 


TWIS program mas OESIG.NED TO TEST FOR SIMILAR CEGREES OF AGREEMENT 
eETwEE'M TmO or more SQUARE ERROR *»*ATRICES 


“I 

NR 


X(I,J] 


a THE NUM0FR OF tables OR MATRICES TO 0€ COMPARED 
a number of rows? alSO TmE NUM0£R OF CCLUMNS SIf4CE THE 
MATRIX IS SQUARE 

a TmE VALUE IN THE MATRIX FOR RQw I ANC COLUMN J 


C 

C 


REAL KHAT,LCL 

DIMENSION x(ao.ao) .sxR(ao),3xc(ao) 

DIMENSION UCL (30) ,LCL (30 ) ,KHA T (30 ) 

LaaO 

MaO 

Kal 

REA0(S,10) ME 
10 FORMAT(ia) 

100 00 aOO lal.L 
SXR(I)aO.O 
SXCf 1)80,0 
00 300 Jsi ,L 
3U0 x(I,J)aO,0 

300 continue 

REAOlS,aO) NR 
20 F0RMAT(I2) 

00 ttOO tsi,NR 

ReA0(5,30) (X (I, J) , Jal ,NK) 

30 FnRMAT(l2(Fa,0) ) 
ttOO CONTINUE 


REA0(S»3l) 

31 FORMAT! 'AAAAAAAAAA AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA A * 


999 


33 


wWITE(6,999) 
FORMATCM') 
wHITE(a,3l) 
mRITE (b, 33) 
FORMATdX, 


' /// 


ARlTE(b, 3 a) 

3a FORMAT!////, IX, 'THE 


35 


wR 




6,35) 


original ERROR MATRIX I5:d 


FORMAT! 1 X, *. 


.*/) 
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II 


ia 

3S 


I? 

IS 

ao 

a 

4 

4 

44 

45 
4 b 
47 
44 
49 

I’’ 

Vi 

5i 

|4 

55 

54 

57 

54 

59 

40 

61 


hi 

63 

64 
45 
66 
67 
64 
69 
7 ft 

71 

72 

73 

74 

75 
74 
77 
74 
79 

40 

41 

Si 

44 

45 
86 


4 io 


up. 450 


coNrxMuc 


XNaO.O 

00 500 t«l»NR 
00 600 jsl*NR 

|xc|j!sixc[j!«x(i,j 

600 CONTINUE 

XNaXN»SXN(I) 


♦XU#Jj 


500 continue 
TH lsO.O 


THato.o 
TMjaO.O 
74430 ,u 
00 700 I31,NN 


rHisTHi*x(i, n 

TH 2 sTrt 2 ^SXR(n * 3 X 0 ( 1 ) 


TH33TH3^X(I, n*(SXff(I)*SXC(I) ) 

00 800 jsi.Nk 

TH43TH4*X(I, J) •(SxR(I)4.SXC(J) J •♦2 
800 CONTIiMUE 
700 CONTINUE 
THlsTHl/XN 
TMasFHa/ (XN**2) 




C 

C 

C 


TH33TH3/ (XN**ai 
TH43TH4/ (XN*»5) 

HMAT(K)3(THl-THai/(l,-THa) ] 

sossQ«r(THi*(i,-tHi)/((i,-TH2)**a)>(2,*((i.-THi)*(a.»rMi*THa-rH3)i 
i)/{ (i.»TMa)** 3 j^(i.-rMi)**a •(rH 4 « 4 ,*THa»*a)/(i,-TMa)** 4 )/xN 


A)/{(1 

the steps that follow calculate the 95X CONFIOE^CE INTERVAL FOR KHAT 

UCL(K)8KHAT(X)-M ,96*S0 
LCL(X)S4HAT(KJ«1 ,96«»S0 
arITE (6,40) 

40 forma T(///, 1 X, *lOwER limit*, 4X, •khAT*,4X, ‘UPPER LIMIT') 

.'•RITE (6, 45) .... 

45 FORMATflX.* , *.ttx.* *.ax.* 

wRITE(6,50) T ( K ) • uCCT^ ) 

50 F 0 RMAT( 3 X,F 8 , 5 , 3 X,F 8 , 5 , 3 X,Fa. 5 ) 

K 3 K *1 
MSM^l 

IF(M.lT.mE) fiU TO lOO 
wRIT£(6,900) 
format (»l *; 


'/) 


900 


wRIT£(6,9lft) 
910 FORMAtax,'- 


SUMMARY table ANO CUMPARISONS ' ) 

•*•***•***•*•***«*•*•*••*'////) 


wRITE(4;92Q) 

920 FORMATdx, ' matrix', ax, 'lower LImIT' , 4X, 'AHAT ', ax, 'UPPER LIMIT') 


wRtTE(6,930) 
930 F0RMAT(1*»* 

00 940 K3l,^ 


%ax, ', 




/»4X, 


///) 


WRITE (6, 950) K,LCL(«),amaT(K) ,UCL(K) 
950 F0RMAT(4X, I2,5X,F4.5,3X,Fd,5,3X,F8,5) 
940 continue 

wRITE(6,960) 

460 forma T(////////) 
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N«M£*l 

OU ISOO 

no laoo j«2*M 


X# (J.LE.I) 60 TO laOO 

IP (LCL(J).6T«LCt(!).*M0,LCL(J).LT.UCLa)) 6C TO UOO' 
iPl |^CL(J)2Lr.yC^(l).AN0.UCL(J).6T.LCL(in 6C fQ UOO 

1000 FoiMAl ( iIJ^MATrIx \IS,* IS significantly OlFFEReNT FROM MATRIX 
A12/) 

60 to 1400 


60 TO 1400 

noo MRire( 6 « 1200) i«j 

1200 FO«»*aT(IX, ^matrix 
AX *,12/) 

IttOO CUNTtNUE 
1300 CONTI.mUE 

rtRlTEIh.lSOO) 

1500 FOHmATI*!') 

STOP 

ENO 


IS not SIGNIFICANTLY DIFFERENT FROM MATR: 


//OATA 
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Appendix IV. Listing of APL Computer Program CONTABLE 


<7COHTA»t.l.HOW CQ39 
9 CONTA»Uh>HOM 

Cl 3 'CONTAeu** 

C23 'ANAC't'SCS or MUI.TXDXMKNS zonal CONTZNOCNCr TA»LCe* 

C33 K , LKK - - - OCrAKTMKNT Or STATS tTXCS» VAX SU • 

C43 •KNTeASP 7/5/1976' 

C53 • • 

C63 ' THIS rrooRAM mxll rsRroAM analtscs or comalkts or zmcomalets' 

C73 ' MULTZDXMCNSXOMAL COMTXMOeMCY TASLCS USXNG A UOOLXNKAR MODEL • 

C33 'AAP-ftOACH,' 

C93 • • 

• DATA SHOULD DC ARRANGED XMTO A COMTINGCMCf TADLE^ AMD MAT DC > 

C113 'STORED IN AN ARRAY SMR RRXOR TO RROGRAM EXECUTION OR MAY DC ' 

C123 'ENTERED URON REOUCST . WHEN AMALY2XHG XNCOMRLETC TADLES THE XNXTXAL' 

C133 'FITTING TADLC MAY DE STORED SIMILARLY IN AN ARRAY NAMED ONES FRZOR' 

C143 'TO RROGRAM EXECUTION. URON REQUEST, THE USER SHOULD ENTER A LOGLZNCAR' 
CIS! 'MODEL WITH WHICH HE INTENDS TO FIT THE DATA. THE LOGLZNEAR MODEL' 

C163 'SHOULD DE ENTERED DY THE CONF ZGURATZONS AMD THE DIGIT 0 USED TO* 

C17] 'SEPARATE CONFIGURATIONS. FOR EXAMRLC , THE NO-THREE-FACTOR-XNTCRACTION • 
C183 'MODEL OF A 3>DZMENSZONAL TADLC, <C12f ^13, C23 ) f IS CMTERCD AS • • 

C193 • • 

C203 • 12013023' 

C213 ’ • 

C223 'THE PROGRAM FOLLOWS AN ZTERATZVC RRORORTIONAL FITTING SCHEME TO ' 

C231 'COMRUTC THE MAXIMUM LIKELIHOOD ESTIMATES OF THE EXPECTED CELL VALUES.' 

C243 'THEN THE PROGRAM COMPUTES THREE GOODNE3S-OF-F IT CHI-SQUARE STATISTICS;' 
C2S3 'PEARSON CHI-SQUARE, LIKELIHOOD RATIO, AND FRCEMAN-TUKEY CHI-SQUARE;' 

C263 'HOWEVER, THE DEOPCE OF FREEDOM ASSOCIATED WITH THESE STATISTICS SHOULT' ' 

C273 ' SS DETERMINED Dt THE USER. AFTER ALL THE COMPUTATIONS RELATING TO THE' 

C283 'CURRENT MODEL ARE PERFORMED, THE PROGRAM ASKS WHETHER MORE HYPOTHESIS' 

C293 '(I.e,, ANOTHER LOGLINEAR MODEL) IS TO DE FITTED, IF NOT, PROGRAM' 

C303 'EXECUTION IS TERMINATED,' 

C313 ■ • 

C323 * all VARIADLCS are localized except SMF and ones, PROGRAM' 

C33;] 'EXECUTION requires the functions named sum AMD YSSNO, all of THESE' 
C341 'PROGRAMS PLUS THIS DOCUMENTATION ARC GROUPED UNDER THE NAME CONTADLEGRP , • 

C353 • • 

C36-J • FOR MOR:E DETAILED DISCUSSION OF THE LOGLINEAR MODEL APPROACH SEE' 

C 373 'THE FOLLOWING REFERENCE;' 

C38] • ■ 

C39] ' DISHOP, r.M,M,, FZCNDER6, 8,E, AND HOLLAND, P,W,, DISCRETE' 

^403 ‘ MULTIVARIATE ANALYSIS* THEORY AND PRACTICE, CAMDRZDOC,' 

C413 MASS,; THE MIT PRESS, 1975.' 

C423 • • 

C431 SAMPLE BUN OF THE PROGRAM CAN DE ODTAINED FROM 3, K, LEE,' 

C443 • • 


*7 
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C13 

C23 

C33 

C43 

Z52 

C73 

C83 

C93 

C103 

CU3 

C123 

C133 

C143 

C1S3 

C163 

C173 

C183 

C193 

C203 

C213 

C223 

C233 

C243 

C2S3 

C263 

C273 

C283 

C293 

C303 

C313 

C323 

C333 

C343 

C353 

C363 

C373 

C333 

C393 

C403 

C413 

C42 3 

C433 

C443 

C453 

C463 

C473 

C483 

C4-J3 

C503 

C5U 

C523 


7C0MTA»LC £03*9 

7 CONTA»LS f D XM } FOLZ> }A|W|T)XM|V)VV| STCA« ) XMP | M ) U ) D ) C ) MNX I MD 

I I 

•THC PCFAULT AAAAMKTBR VALUKt rOA TH8 XTCRATXVK AAOAOATXOMAL ' 
•rxTTZNo scHKMC aak; MAKXMUM PKVXATXON m 0*01* 

' MAXIMUM NUMPCA OF XTCRATZONS m IS*' 

I I 

»C0XM:>KNTKR THC DXMCHSXONS OF THC TA»CB * 

OXMf-a 

•ZS THC DATA ALACAOT ZN AARAT tMFT(YCSf NO* OA CTOA ) • 

W4.TCSH0( A^Q) 

•♦( (W«0) f (W«l) )/0»»o 

ACSTAAT: • CHTCA THC DATA WZTH LAST SUDSCAZAT CHAHOZMO FASTCST • 
SMF«.Q 

eCTt'AMl MOAC DATAT' 

W<- r eSHO( AlQ) 

-»( (WaO) » <w«2) )/0»oo 
'ENTER MORE DATA* 

3MF4-SMF, T«-Q 
•*6ET 

GO:^( ( x//»SMF)sX/DZM)/001 

•ERROR. HUMSER OF ELEMENTS ZNCORRCCT. HCAB IS THC ZMRUT • 

SMF 

•DID rOU ENTER ZERO COUNTS? LETS TRY AOAZN' 

■4RCSTART 
OOl :SMF«-DIMf SMF 

■ DO YOU WANT TO SEC THC TABLET • 

W<. YCSNO( A«.Q) 

■»( ( WaO) f (W«2) )/0»OOO« 

SMF 

goon: 'ANY FZXCD ZEROS?' 

W*- r ESNO( A<-C1) 

■♦< (W»0) * (W«2) ) /Of soil 

•ZS THE MATRIX OF ONES AND ZEROS ALREADY IN ARRAY ONES?' 

W4. YESN0( Af-Q) 

-»( ( w*0) f <w»l ) )/0f 007 

•ENTER A MATRIX OF ONES WITH ZEROS IN THE ARRROFRZATE FLACES > 
G022:ONESf-a 

006: 'ANY more?* 

W«. YCSNO( A«.g) 

-»( (WaO) f (W»2) )/0fOO7 

•ENTER MORE ONES AND ZEROS' 

ONBS<-ONES*a 

•*006 

G 07 t A4.f « ones 

•»< ( < X/OIM)*A) , ( ( x/DIM)<A) ) /OOl 9 *0020 

• tOU HAVE NOT ENTERED ENOUGH NUMBERS rCT' ALL THE REMAINING* 
•NUMBERS ARC TAKEN TO BE ONE* 

ONESfONCS* V ( x/DZM)-A)f 1 
40021 

0020:‘'OU HAVE ENTERED TOO MAN t NUMBERS;* 

0021 tONCSADZMf ONES 

•HERE IS THC MATRIX OF ONES AND ZEROS AS IT STANDS NOW 

ONES 
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C333 

•ZS THIS TMK MATKXM rOU WAMTT • 


CS4J 

W4.'rKtNO( A4>0) 


CSS3 

4( (W«0) > (W«l ) >/0»»O2 


C563 

•RK-KNTKR TOUA MATRIX OF OMCt AMO ZEROS* 


i:S73 

40022 


C583 

00l?:0NB8^DZMf ONES 


C593 

•DO TOU WISH TO SEE THE TADCS OF ZEROS AND 

ONES?* 

C603 

W«.rBSNO(A4.a) 


C613 

4<(w«0)»(f*2))/0»oo2 


C623 

ONES 


C633 

4002 


C643 

OOll :0NES4>DZMf 1 


C6S3 

002: 'ENTER MAROZNAUS TO >E FZT* 


C663 

ooi 8 :v«-o 


C673 

CrCLEt'ANT MORE MAROZMAI.ST' 


C683 

W4- t SSNO( A«.Q ) 


C693 

4( (W«0) » (W«2) )/0»OO3 


C703 

■ ENTER MORE MAROZNACS* 


C713 

V*.V,n 


C723 

4 CTCI.B 


C733 

003 : v<.v,o 


C743 

I 4-0 


C7S3 

J 4 -O 


C763 

GO 12 : 


C773 

4< J>#‘V)/0O13 


C783 

4(vcj3«0)/«o12 


C793 

vv^,vcj3 


C303 

0014:<J4>'J4l 


C313 

J>fV)/0016 


C823 

4(vcj3*0>/ooi7 

•“igj? 

C833 

vv«-vv,vc J3 

C343 

400J4 

C353 

CS63 

©Oi7:4(VCj-i3.0)/OO14 

I4>I41 

C873 

A«.p VV 


C383 

VV«.VV, < < I<fCiXM)-A)/-0 


C39] 

40014 


C903 

GOl3:*rOU HAVE MOT CNTCgED AH’t MAftOXNAUS T CT , TRY AOAZM * 

C913 

•»ooi 8 


C933 

GO 1 6 : ( I f f D I M ) f VV 


C933 

(v'r'OU HAVE ENTERED <),<^Z),^' MAROZNACS.' 


C943 

HDf-O.Ol 


C953 

MN Z ». 1 5 


C963 

FOCD4-ONCS 


C97] 

STBFS<-I 


C993 

M«-l 


C993 

* DO tOU MISM TO SET THE ^AAAMETEH: VALUES rOA THE XTSPATXVE 

Cl003 A«-Q) 

C1013 -*( ( w«0) » (w«2) )/0»»xo‘-oo^ 


C102 

3 005 : 'WHAT IS rouR vacue of maximum DEVXATIONT* 

C1033 

Cl 04 3 'WHAT It rOUDt vAuuc or ma:<zmum NUM»CP: of 
C1053 

tTEF:ATIOW«7 * 

Cl 06 

• S I04.00R ; ^ < M > MNI > /NOCON 
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C107] c<-0 

C106] t«-l 

C1093 UXTLOO^:-*( t>«TBI^«)/CHKCK 
CU03 XN&<.(*(vcti3«0) )/vCXO 

C1113 ^OL.D«>rOi.X»X ( XMD «UM tMr)*<(XNO sum FOCD)-f(XM» «UM FO(.O)|0> 

C1123 

C1133 -*uXTuoof- 

C1143 CH*c«:XM 

C11S3 CHCCNCH!:4( X >STC^S)/OUT 

C1143 *M»«*<*(VCXI3«0) )/vCX»3 

C1173 C 4 >< r/ I » ( XNX> SUM SMr)«(XMD SUM rOCO ) ) 

CU83 -»<c<o>/ooio 
CU93 ix-c 

C1203 Ooio;-»(C>MD)/STB^ 

C1213 

C1223^ ^CMCCicBft 
CI 233 ' »TKPJI4*.M>1 
C1243 ■^fiai.oop 

C1233 ooT;u«t( ( ,FOU»gO)^( ) 

■■ 

C1283 (^'MAXIMUM DCVXATXOM a ' ) t ( ) » ( T ' MO, OF XTHftATIOM a •)»♦*« 

C1293 • • 

C1303 • • 

C1313 (♦•FCAftSOH CHX-SOUAFC • ) ,^^/( ( ( •U)/( ,SMF-FOUO) ) •2)-( 

C1323 r<-<SMF*o.5) + ( (»*»F+i) ,o.S)-< (l+4xFoi.i*)t0.5) 

C1333 (♦•FFKKMAN-TUKKf CHI^SOUAFC ’ ) » f»-/ » ( T« 2 ) 

C1343 u«.,SMF,o 

C1353 (♦•CIKKI-XMOOD FATtO < ) , ^2* ♦ / ( < ) / » *^F ) x • ( ( »U ) / , SMF ) - ( ( /wU ) / , FOUf ) 

C1343 (♦•totac mumscr of Fxxeo scfos • > x/x>xm).<.«./»ohcs) 

C1373 * • 

C1383 • • 

Cl 393 'CO 'OO MAHT TO SEE THE C;<FeCTEI' TAS4.CT • 

Cl 40 3 W«-«BSHO(A«.fl) 

C1413 +( (»» 0 ) » («« 2 ) )/ 0 »coNT 

C1423 ■ ■ 

C1433 • • 

C1443 ( L0.5^«'o»-cxi00)^100 

C1453 • • 

C1463 • • 

C1473 ‘CO rOU WAHT TO SEE THE FFEEMAH-TUKE'f S>BV X AT XOHST ' 

C1483 W4-<ESH0(A^Q) 

C1493 <w» 0 ) » ( w* 2 ) )/0»oo9 

C1503 • • 

C1513 • ’ 

C1323 (LO.S^T*lOO)^lOO 
C1533 • • 

CIS43 • • 

C1SS3 Q'^Ot’DO rou WANT to see the FESIDUACS <06SEFVED - S::FSCTSt> > 7 ‘ 

C1563 W*- rESNO( A«-0) 

CIS73 •♦l (M»0) » (M»2) >/0 .go4 
C1583 • * 

C 1593 * • 

C1&03 (L0.3t<»mf.foi_&)x100)-^100 

C1613 • ' 
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CU23 ' ' 

C163] 004; 'C*0 WANT TO SCC THC «rANOAA»tSKO ACSZOUAkST' 

C1A43 W<.rK#NO(Ai.0> 

C16S3 •»< <w*0) , (W.2) )/O»®O0 

C1663 ' ' 

CU73 ' ' 

C1683 (L0.5+< (•MF-F0U0)*( (FOuo+roLO|0)*0.S) )«100)T100 

C1693 ’ • 

C1703 ' ' 

C1713 009 : 'DO rou WANT TO «CK THC LOO C»ACCTATZONS? • 

C1723 w^retNO(A^0) 

C1733 ■»< (w«o) f (w« 2 ) )/ 0 »coNT 

C1743 • • 

C1753 ' ' 

C1763 (L0.5T<o^OLD+roL»i0)xl00)+l00 

C1773 ' ' 

C1783 • • 

C1793 CONTt'WANT TO Tftf SOMC CTNCR MfPOTMOKST* 

C1303 W<.l ClHO( A*.Q) 

C1813 -♦(< waO' »(“■!) )/ 0 »oo 2 

C1323 ’ WANT TO ANALTZE ANOTHER TA»LET ' 

C1833 w«-fE»NO( A«. 0 ) 

C1843 -♦( (««0) » (W*! ) )/0»»osiN 

Cl 83 3 

C18o3 N0C0N:(^'N0 convergence ARTER •),(^HNX)t^< ZTERATXONS^* 
C1373 (^'TME DEVXATIOfJ AT TWX» OTAOC X» * ) # ♦» 

C1883 'CO rou want to do more ITBRATXON* or change iOUR VAt-UC* 
C1893 manxhuh oeviatxont' 

C1903 w^rE*NO(A«.Q) 

C1913 <w« 0 ) , (W.i ) )/0»oo5 

C1923 ■♦CONT 
C1933 SHC; ' • 

C1943 ' ' 

C1933 •• 

T 


VCONTINOCHCl CD3^ 

9 CONTINGENCt 
Cl 3 • ENTER OBi«ERVEr> CELL RREOUENC lES »f ROWS' 

C23 Oi-\0 

L33 h;o4-o,q 

C 43 -«U 1 t I ' T • «( n,Q«- • more OSSERVEI- RRCOUENCIES ro CNrCR-T ' ) Cl J 

C5J 'ENTER NUMSER OR ROWS ANE' C0LUMNS« RESRECT I VEL > ' 

CA3 DR^X/<SHARE*.a)-l 

C73 ROWTOT,-4./MAT^SWARBrO 
eg] IH.^/COuTorr^^MAT 
C9] EUR ♦.ROWTOT. , kCOUTOT-N 
C103 CHISOUARE*.*/ , ( (MAT-EKR ) t2)TEXR 
CllJ ■06'SERVEC' rBEOUENCIES ARE; ' 

Cl23 

C133 'EMRECTEt' RRBOUENCIES ARE; 

C143 e'tP 

ClSl ' » ' •-•'ONTINOENC t CHX-tOUARE WITH • ) , ( ' 

9 

\ 


C'R r i 


. • ♦CHI'iCfUARE 


48 


9TWOWAT CQ37 

9 O TWOMAr r;A|C; K )M)R}S)SXDE)T)TAC>|TOP)TV}V 
C13 TOP^Q^OfQl* ' top HEADXMO' 

C2] SIOBf-QrO'^QT ' enter SIDE HEADING* 

C33 T4.DC23+Oxst-oci3+OxR«“CC33+OxC«.DC43 

C43 TAB«.(S,T)fOXN«-l 

C53 

C63 v<-( fC»R]«H)/\590 

C73 T2:TAeCM>K]f.+/fCV}C3aK 

C8] -»T2x »T>.K«.K + t 

C9] -kTlx 
C103 " 

cm '»Top 

C123 

C133 94-PTV4-* •.<9xT)f*110 

C143 TV FMTtT 

C153 

C163 SIDE 
C173 TV4.84.TV 

C183 «T1 

C191 t3: • I • 

C20] (T* •)»(T«)»(T' I ■ ) ,^TV FMT TABCKn 

C213 -»t3x »s>Kf.K + i 

C223 

C233 QTa4-0<-'' 

C24] (t'FOF. THIS TABLE. TOTAL RESPONSES APE; • ) , ^+/ . TAB 
O 

\ 


otesno CD39 
9 W 4 >rESNO A 
Cl] A4-4rA,* 

C23 ■♦(( A/Aa ' STOP •),( a/A= • YES • ) , ( A/Aa • NO •))/Gl»02»G3 

C3] ‘TYPE YES, NO, OP STOP PLEASE* 

C41 A,-4f ( A«-a) , • 

C53 T2 
Z61 oi;w4-o 
C7] 40 

C33 g2:w«.i 

C93 40 

C103 g3*w<-2 

9 

\ 
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7TMRKi:WAr CQ39 
7 W THASBWAY ;; 

C13 DEp-TM^-OtafOtO** ' DEPTH HEADXHO* 

C23 Kow^-stOf OtQT • gnter row hsadxmg* 

C3:i COi.UMN«.Qr0tQT ' CNTEP COLUMN HSADXHO* 

C43 o<-wci3>0*'*=TWC23+0xCt-w(;334>0x5<.wC43T0xT*.wc53T05«'J<-“C63 

CS3 TAD^(O,P,C)f0XtTlT0X««-lffX 

C63 Ti ; *4-1 tfY«.MC(«C»* 3«x >/»<»» 3+0xJ«-l 

C73 T2:2»-rc{YC»T].j)/\s»3+0xK«.i 

C83 T3:tadcxj j>K3«.4-/2CI03aX 

C93 -»t3x iCiK^-K + i 

C103 -»t2x \ ft > J4-J+1 

LllJ -^rix 

C;i23 OtOt' • 

i;i33 • •, COLUMN 

C14] •• 

C153 s,-fTV 4 .« ‘»<9xC)fiiO 

C163 TV fmt\C4.0x1«-1 

C173 

C183 DEPTH, ROW 
C193 (ff)f-* 

C203 TVM34TV 
C213 T4:j«-1 

C223 t3: • 1 • 

C233 (7' 111 111 I' PMT( I, J) ) ,»TV PMT TABCItJfl 

C243 tT5x \ Pi'JTJ+l 
C2S3 7T4X > Diit-x + i 
C2A3 

C273 Qt-a^QR'* 

C283 < ♦ ' POP THIS TADLE, 

7 


TOTAL RESPONSES ARE 


> , 7+ / , T AD 


so 


Appendix V. 

A Quanclcacivtt Machod co last for Similaricy 
bacwaan Photo Intarpratars 

Russell G. Coogalton, Research Assistant 
and 

Roy A. Mead, Assistant Professor 
School of Forestry and Wildlife Resources 
Virginia Polytechnic Institute and State University 
Blacksburg, Virginia 24061 

Biographical Sketch 

Mr. Congalcon is a native of Mew Jersey and received his B.S. 
in Forestry at Rutgers University. He is presently a graduate research 
assistant at V?I and SU. Dr. Mead received his B.S. in Botany at 
Northern Arizona State University, his M.S. in Reinote Sensing at Colorado 
State University, and his Ph.D. in Remote Sensing at the University of 
Minnesota. 


Abstract 

A method has been developed to quantitatively test the degree 
of similarity between photo interpreters. This method involves giving 
each photo interpreter the same set of photos to interpret. An error 
matrix is then . generated for each interpreter by comparing his interpre- 
tation to the actual ground cover. This error matrix is then analyzed 
using a computer program called KAPPA. This program uses discrete multi- 
variate analysis techniques to determine if one error matrix (i.e. , photo 
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laeerpracsr) Is slgnificsntl/ dl£f«r«nc from anothmr. Ihs program 
can ba alcared co taac for slmllaricy ae dlffaranc confldanca lavals. 

Not only doas this tachnlqua allow ona to compara two saparata in* 
carpracars, but it also allows ona to cast vbathar an individual pboCo 
intarpratar is eonsistant through tlma. 

Introduction 

Photo intarpratacion is the art and science of identifying 
objects and deducing chair significance on aerial photos. Good, con* 
sistant photo intarpratacion depends upon the experience and skill of 
Che individual who dalineatas the boundaries becvaen vagetatlon/land 
cover types over the landscape. The judgment involved is generally 
qualitative in nature, and therefore difficult to evaluate or compare 
with Interprecatlons made by others. Usually the interpreter has 
intuitive feelings about how well he is doing, but is unable co support 
these feelings with any specific tests. This paper suggests a way of 
quantifying photo interpretation results and g±res a statlsclcal method 
for comparing these results. 

The procedure proposed in this paper can test for Che degree of 
similarity between interpreters, or test the consistency of the same 
interpreter over time. Testing co see if interpreters are similar is 
useful when more chan one interpreter is to work on a project. If it can 
be determined that the delineations made by all interpreters are not 
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significantly dlffarent, chan cha project will yield uniform results 
for all Interpreters. Also, It would be useful to test the same 
Interpreter over a period of time to check for changes In his Inter- 
pretation. It nay also be important to determine If varying types of 
photography (fllm/fUcer combinations), or seasons of photography result 
in significantly different delineations. By placing a grid over each 
delineation, the individual cells are assigned to the land cover/vegecaclon 
type which represents the majority of the cell. Each cell is then com- 
pared one-by-one with the corresponding cell (l.e. , in the same location) 
from another delineation. If one of the delineations Is assumed to be 
correct (reference data) , then comparison of the two secs of spatially 
defined cells yields a measure of "photoincerprecation accuracy". This 
is usually expressed in the form of an error matrix. 

Procedure 

An error matrix Is a square array of numbers sec out in rows and 
columns which expresses the number of cells assigned by the photo interpreter 
to a particular land cover type relative to the actual land cover (reference 
data) . The columns represent the reference data and the rows indicate 
Che photo interpreter assigned land cover type (Figure 1). 

Reference Data 
ABC 

Photo ^ 

Interprec- 
caclon 

C 


! I 

I I 

T i 


Figure 1. Error matrix format for three land cover types. 
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Th« Qu&b«rs la ch« arror macrtx are telllee compiled by com* 
paring the photo Incerpreeatloa with the actual cover typ^ (reference 
data) on a cell by cell baala. All correct clasalflcatlona are 
located on the major diagonal of the error matrix. 

The specific method used to generate an error matrix Is dependent 
on what Information Is needed. If the degree of similarity between 
two or more photo Interpreters Is to be determined, each Interpreter 
Is given the same aerial photographs to Interpret. An error matrix 
Is then tabulated for each Interpreter by comparing his Interpretation 
with a reference data sec (correct delineation). If the test Involves 
determining the consistency over time for a single Interpreter, Chen a 
representative part of a selected stereo pair Is Interpreted at the 
beginning of a project. Ac some later dace the remainder of the photos 
are interpreted and then the two error matrices (Time A and Time B) 
are compared. Finally, If it Is desired to measure the accuracy of 
delineations made on different types of photography, a separate inter- 
pretation is performed on the same area for each sec of photos by each 
interpreter and an error matrix Is generated. 

Once Che error matrices are generated, a discrete multivariate 
analysis procedure (Bishop ^ al . , 1975) is used to test the degree of 
similarity between the error matrices. This test Is based on a maximum 
likelihood estimate of the multinomial distribution (Equation 1). 
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EX- I (X * X ) 
!■! ^ t«l ^ ^ 




wh«ra: 


( 1 ) 


K •• # of rows In nacrlx 
X^^ * # of obs in row 1 and col. 1 
X^^ ■ narglnal total of row 1 
X^^ ■ marginal total of col. 1 
N • total # of observations 

This equation yields a value RHAT which Is a measure of the actual 
agreement minus the chance agreement. A confidence Interval at a given 
a- level is then placed around the value of RHAT calculated for each 
error matrix. If the confidence Interval for one error matrix overlaps 
the confidence Interval for another error matrix, the two matrices are 
said to be not significantly different at that a-level. However, If no 
overlap of the RHAT confidence intervals occurs, then the matrices are 
said to be significantly different at that s-level. 

This entire comparison process can be performed using a FORTRAN 
computer program called RAPFA. Given the error matrices to be analyzed, 
the program calculates a RHAT value and a confidence interval for each 
error matrix. The program then prints out which error matrices are 
significantly different and which are not. 
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Th« data usad In thla study wars takaa froo Lauar at 
(1970). Flva photo Intarpratars intarpratad tha sama aarlal photo- 
graphs of Yosamlta Vallay, California, and thalr Individual arror 
matrices ware ganaratad. Also, flva flla and filter combinations 
ware usad with a single interpreter, and arror matrices ware generated. 

Results and Discussion 

All five of the Interpreters tested on the photos from Tosemlte 
Valley produced significantly different delineations (Table I), the 
confidence interval was calcxilated at the 95Z level. 


Table 1. Summary table for five Interpreters of Yosemlte Valley 


photos. 

Interpreter 

Lower Limit 

KHAT 

Upper Limit 

1 

0.31167 

0.31991 

0.32815 

2 

0.28623 

0.29420 

0.30216 

3 

0.36677 

0.37485 

0.38293 

4 

0.23115 

0.24156 

0.25197 

5 

0.20878 

0.21925 

0.22972 


The results of the five different film and filter combinations 
are presented In Table 2. These results were also calculated at the 


95S 


confidence level. 
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Tabl* 2. Sunnary tabla for eha flva flla and filtar comblnatlona. 


Fila/Filcer 

Lower Limit 

KHAT 

Upper Limit 

m-301/W2S 

0.3U67 

0.31991 

0.32815 

IR/W89B 

0.29615 

0.30436 

0.31258 

Beta Aero IR 

0.11318 

0.12071 

0.12825 

Enhancement 1 

0.25427 

0.26163 

0. 26898 

Enhancement T 

0.36704 

0.37438 

0.38173 


AS can ba sean from Tabla 2« cha Interval for IR*301/W25 over* 
laps wlch the interval for IR/W89B. Therefore, these two interpre* 
cations are not significantly different. All the ocher intarprecatlons 
are significantly different. 


Summary 

The examples given in this paper indicate how photo incerpre*' 
cation results can be quantified using error matrices. These error 
matrices can then be compared using a discrete multivariate analysis 
procedure and conclusions made. 
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‘ ABSTRACT 

A vorkiat ceaftrmet was ia Sioux Folio, Soutk Ooketo Novoabor 12, 
13, oad 14, 19S0 dtoliat uitk loatfsot Clotoifieotlea Aeearoqr Aooossaoat 
Ptoetdurto. Tkistooa foaol proaoatocioaa utrt aoda oa throo loatrol topics t 
(1) ooapliat procaouros, (2) stotistieol oaolytis tockaiquos, oad (3) oxoaploo 
of pcojocts vkiek iaeludod oeeuroep oososssoat oad tho ossoeiotod coots, 
lotiotieol prebloao oad ooluo of tha occarocr dots to tho roaeto soaoiaf 
spteiolist oad tho rosoucco aoaottc. Noocly twoaty eoaforeact ottoadoos 


pocticipotod la two discusoloa sosoioa oddctsoiat ▼orious lotuts ossoeiotod 
with oeeuraey ossossatat. This popos presoats oa oceouat of tho oceonplish< 
aoats of tho coaforoaeo. 
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nmnDucnoM 

Xa tha Ttara siaet laadiat iMfarr first bseaat svsilsbls» «a uatold 
aiabar of Uadsot tesaas turo boaa di|i tolly oaolysad to elossify load eovas. 
Tbata cloooifieotioaa oca oat vithoat actor, oad boaa baaa siabjaet to elooar 
sccatiay by erities oad potaatlol oaara thoa aiallor prodaeto daaalopad by 
aero tcoditioaol aathods. A fav potaatlol ooars of Loadaot doto wara diocoaragad 
by tha iiafalfillad oxpaetotioaa splcltad by tba rasalto of aorly iavasti|o- 
tloas. This bos raeaatly led resasrebars sod goaamaaat agaacias to proceed 
eotttiouoly with tecbnology traasfer. Thus, scieatists bova bean keenly swore 
of tbe need to sssess tbe aceucscy of Loadsot elossificotioas before dis* 
tribatiag tbe products to users. 

Topegrspbie aoppiag procedures iaeluda routlaa aaoluotioas for co^plioaca 
with wall defiaad seeuraey stoodords sad tba ocearsey sttoiaobla tadar specific 
coaditioas (tarraia cborscteristics sad aoppiag aquipSMnt used) ora %fall 
kaowa. This capability Is tba result of aaay directed resaorcb efforts. 

However, teebaiquas for sssessiag tbe scearocy of loadsot elossificotleas have 
developed la sa od hoc aoaaer. Huef such aatbods era aot statistically souad 
sad can yield biased estlsutes of aceucscy. 

For exaople, researcbers used tba llaited available ground lafonation 
(l.e., asps, pboto lacerpretatlons or less often actual visits to tbe field) 
collected for developmaat of traiaing statistics to estloate classlflcatloo 
accuracy. This can result In over optiaistlc estlaates of classification 
perforoaace, particularly when tbe training data does net adequately describe 
tba scene variability. Wlndsbieid surveys, la which a few easily accessible 
arass are visited oa tba grouad, are oaotbar biased approach to accuracy 
assessaeat. In addition, biases can also be Introduced by ualag s different 
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classif ie«tiea fr a— work foc» aceuraqr Ms«saMat thM tkac as«6 ia 4t— loplag 
tkt di|iul clasalfieatiea. 

Tha krtad, aert raetatly, baa baaa to aaapla tha elaaaifieatiou aad aaaaura 
tka dagraa of afraaaaat vitk a aat of apacially dafiaad rafaraaea data (i.o., 

grouad truth). Aaalyaia of tha raaultiag accuracy data can guida raaaarehara 

0 

ia acaaa salactioa (aaaaoa, ate.), and datacaiaa tha aoat appreprlata —thoda 
of elaaaifleatioa for particular appUeatioas. 

Tha iaportaaca of aaaaaaiag elaaaifieatioa accuracy, tha lack of aay 
ataadard procadurta, aad tha liaitad aiiabar of raporta ia tachaically rtvitvad 
Jeuraala, Juatifiad tha eoaftraaea diacuaaad ia thia papar. Oely a ralatlvaly 
aull auabar of raaaarehara haaa workad ia tha aubjaet arta to aay great 
artaat. Tharafora, attaadaaea at tha eoafaraaea waa liaitad aad by iarltatioa 
oaly. Tha apacifie objeetiaaa of tha eoafaraaea ware: 

1. To datezaiaa tha atata^f*tha*art for accuracy aaaaaaaeat procedurea. 

2. To provide a forua for exchaaga of idaaa coacaraiag accuracy aaaeaa* 
aeat proeadurea. 

3. To idaatify reaaarch aaada aad racoaaaad tha approach that ahould 
ba taken to iaprova accuracy aaaasaaaat procedurea. 

C 0 1 TOE3 I CE TBEHES 

A coapreheaeive proceediaga of the fomal eoafereace preaaatatiooa ia 
planned. How«v«r, it ia worthwhile to identify and auanarire the aajer theaea 
that dc-veloped froa the conference in general. 

Accuracy ia a aaaaure of the aaouat of agreeaent between two data seta. 
Typically thia is a theaatic aap ia gueatioa and a reference data set often 
thought of as ’’ground truth." Bowavar, whan this procedure ia generalized 
other applications becoaa apparent, including change detection analysis for 
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tlM Maltoriag of particolar rttourett. Furthocaoro, toquoatiol approisol of 
a elasslfeatioa caa ctsalt ia battor oad cMules. 

Thera are saaaral type* of aeearaey aad it is iaportaat to idaatify ahich 
ia baiag utilitad. Two aajor eateforlas of aceoraqr are aita specific aad aoa 

aita apacific accuracy. Noa aita apacific accuracy coaparaa tabular auaonrya 

* 

of the proportioaa of the area aappad iato each of the catagoriaa. Site apacific 
accuracy utiliraa the apatial aatara of the data. That ia, two apatially 
defined data seta are regiatared and coapared for the aaouat of agreeoaat. 

This can be a polygon, grid cell, or point cooparisoa. In this caae, the 
difference between the two data sets results in a spatially defined binary 
data set. This represents the population we are saoipling for the paraoeters 
in question. 

An error aatriz or contingency analysis approach to accuracy assessment 
is still another method of comparison of the two data sets. This requires a 
site specific (spatially defined) approach. 

Fttrthesaore, many factors affect the validity of an accuracy assessment. 

The quantity aad quality of ground truth depend upon the aethods used for 
sample size determination and data acqusition. In light of this, it becomes 
apparent that the term **ground truth” is ill defined. What is "ground truth" 
with regard to parameters such as percent of ground cover? Caa this ever 
really be measured? For many cover types, this parameter can be estimated 
more accurately on aerial photographs than by ground procedures. 

Finally, one should not lose track of the difference between the use- 
fulness of a specific product and its estimated accuracy. A numerical report 
of product accuracy smy say nothing of how much use the product gets or how 
well it compares with what was previously available. A quantitative accuracy 
assessment resuxts in a numerical sumoMry which nay or may not represent the 
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vutfulMtt of tho product. la pooy iaataaeot, « elottifieatioa of low or 
latoraodioto oeeuraep is a waleoat and usaful product. 

Hm dasirtd laforaatioa as wall as ehs aaturt of ths scaaa which was 
clasaifiad* datatuias which is tha aost appropriate oesas of assessiag accuracy. 
Cartaialy, differaat laadscapas aay aaad to ha sauplad diffaraatly for bast 
results. Thacefora, studies should be dona to look at tha saaaitiwity of 
accuracy eatiaatas whaa differaat saapliag procaduras are used. Xa coapariag 
aad assessiag saspliag procedures for accuracy assessaeat, aot only is statistical 
variability to be considered* but also the spatial diversity of the data. 
Furthetaore* all of these considers tions Interact to deteraine the aost appropriate 
saapling and estiaation procedure to use. Much work reaains to be dene* 
utilising designed ezperiaenta with specific hypotheses* to identify tha 
relative reliability of various saapling procedures. 

Assessiag aad reporting* by soae standard oeans, the accuracy of a theaatic 
classification will become aore vital as these products becoae a part of 
geographic inforaation systems. This will be necessary to insure hi^ quality 
output products and well infotaed aanageaent decisions. 

The use of training data for accuracy assessaent results in a soaewhat 
biased but possible useful estimate of overall accuracy. The nature of the 
bias is to overestimate accuracy. The amount of bias depends upon bow well 
the training data represent the variability present la the scene. In soae 
instances, such an approach will be adequate. However, for close scrutiny and 
for within class estimates of accuracy, and independent accuracy assessaent is 
warranted. An approach to aiaiairing the cost of an independent accuracy 
assessaent is to collect accuracy assessaent data at the time the training 
data is collected. This data should be earmarked for later use and not used in 
the training process. 
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AltiMtttb thia eanfartaca did meh to tatoblish eeanmieotloo aaoat rotoareh* 
on utiliaiai oeeurocy ooMooaoat proeodurtti wteh verk remoliu to bo done la 
suoBariziag what procodurea are aoat cooBoaly utlliaod. la additioa, a bibliography 
of the litoratttco aad available cooputer prograaa should be coi^iled aad 
published. A survey of cesearchesa ia the field will help to defiae how well 
they caa aap warioua cower types. This will assist ia dewelopiag a set of 
■appiag staadards. Although accuracy requireeeats aay wary aaoag cower types, 
acceptable aap accuracy staadards are aeeded to OMteh iateaded uses. Staadards 
such as **secoad order at level II** esa help ia aiaiaisiag subjective ewaluatioas 
aad fixully, perhaps aaay classificatioas are sore accurate thaa we thiak due 
to geouetry probleos aad edge pirels. It becoees apparent that elassificatioo 
error aad napping error are not oae ia the sane. Much work seeds to be done 
to discrifliaate between the two sources of error. 
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SlMHMff 

Many iasuM ware diaeuaaed aad debated by the partleipaata. Topics for 
further research were identified and najor theaes sunaarized in this paper. 

The participants recooaended that a working group be established to write 
a '*aanual" or "guide book" on accuracy aasessaent procedures. Possibly this 
group could be focaed as an ad hoc coaaittee within the Aaerican Society of 
Photograaaetzy aad seek funding to prepare the docuaent described above. 

Plans are now being aade to do this. 

The conference succeeded in accoaplishing the three objectives stated 
earlier. A eoaprehensive proceedings is planned tihich will represent state'of* 
the-art accuracy assessment procedures. 
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Appendix VII. 


Listing of FORTRAN Computer Program WILDINT2 
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.A COMPUTERIZED SPATIAL .\NALYSIS SYSTEM 
FOR ASSESSING VTILDLIFE HA3IT.AT 
FROM VEGETATION MAPS* 

Roy A. Mead, Terry L. Sharlk, 
Scephen P. Prisley, aad Joel T. Helnen 
Deparcaenc of Forestr;^ 

Virginia Polytechnic Institute 
and State University 
Blacksburg, VA 21061 


ABSTRACT 

Vegetation and land cover patcems affect the quality of habitat available 
for vrildlife. Given the degree of Intersperslcn of cover types and relative 
value of each edge rrpe and the iaportance of spatial diversity, an inde;: of 
habitat spatial diversity can be computed for each parcel of land (of 
any desired size) relative to each vildlife species or group of species. 

This is accomplished by defining a grid vhich is either placed on a land 
cover 3ap or on an aerial photograph. Each cell is then coded on the basis 
of (its predominant) cover type. A computer program subsequently analytes 
the arrangement of these coded cells and produces maps of (a) interspersion. 

• Presented at I'th .\nnual Meeting American Society Photogrannetr;?’, 'oashingto 
D.C., February 23-25, 1981. 
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(b) juxcoposlcion, and (c) spaclal dlvccsicy. S«paract siulclcolor aaps 
can be aade for any wildlife hablcac of Inceresc using « dlgical fila 
recorder. These asp overlays can be used by the resource manager to 

0 

compare wildlife habitat quality and potential with aaps for forest, range, 
watershed and recreation potential. 

imODUCTION 

There is a tremendous need to develop quantitative aethods to assess 
wildlife habitat. This was specifically aandated by the Resources Planning 
Act, as well as other legislation. Wildlife habitat must be considered in 
all aanageaent plans together with tiaber, range, recreation and watershed. 
While timber inventories have been conducted for aany years, techniques for 
quantifyiag the wildlife habitat still need to be developed. 

The technology of remote sensing has provided the means for aapping land 
cover/vegetation over ver:? large areas for wildlife habitat aanageaent 
(Pengelly, 1973). However, t.he aaps themselves only partially fulfill the 
inventor;? data needed by biologists who must aanage far wildlife. The 
aaps must be analyzed and interpreted to enhance the various characteristics 
of the landscape which have a bearing on aanagement decisions. In short, 
the standard land cover map is a source of information chat aay be helpful 
in aaking aanagement decisions. 

This paper suggests a means to analyze and interpret aaps of land cover 
to produce spatially defined data that will be valuable information for 
managing wildlife habitat. laphasis is on the tec.hnique and not on the 
controversial issue of defining habitat QC.ALIT;. It aust be understood 
that the landscape characteristics iaportant in habitat evaluation vary 
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according to region and the specific wildlife species of interest. The 
various weighting factors discussed in this paper nust be deteminei by 
wildlife managers familiar with local conditions or from agency handbooks 
which give the habitat requirements and preferences for many species. 

The specific objective of this study was to develop a computerized 
system for measuring the spatial diversity component of wildlife habitat 
from vegetation maps . 

Study .Area and Input Data 

The area used as an exaaiple for testing the wildlife habitat analysis 
techniques described in this paper was the Great Dismal Swamp. The area 
is managed by t.he U. S. Fish and Wildlife Ser*/ice as a game refuge and 
includes approximately 34,900 hectares. This wetland was thoroughly 
described by Garrett and Carter (197"’) . The area was ideal for evaluating 
the proposed habitat analysis techniques for three reasons. First, 
the Dismal Swamp "contains a remarkable diversity of vegetative communities" 
(Garrett and Carter, 19'”’) . Second, the area had recently been mapped 
(Gammon and Carter, 1979). Third, the local resource managers were 
available for assistance in evaluating the validity/usefulness of the 
final habitat quality maps that were produced. 

The vegetation maps produced by Gammon and Carter (1979) contained "4j 
separate canopy designations and 143 specific vegetative communities...". 
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This ssp tfss ovtrlaid wich s squirt grid syscin oriincsd in i North-South 
ainnar. Each csll containtd 22 hactaras and foraad a aatrix of 93 rows 
and 42 coluaos. Each call was givan a cotnniunlC7 dasignation according to 
tha covar typa which occupiad tha aost araa within that call. This informa- 
tion was storad on discs for analysis by tha cooputar. 

Given tha above data tha following procedure was used to assess wildlife 
habitat dlvarsit;/' for the Great Dismal Swamp. Since the priaar:^ thrust 
of this paper is to present a proposed technique, all additional inputs 
(e.g., juxtaposition weighting factors and restrictive factors^ are purely 
hypothetical, as is the selected wildlife species "A". 

HABITAT ASSESSMEMT PROCEDURE 

There are four components that fora the package of techniques used for 
assessing wildlife habitat: 

1. Input data 

2. Measurement of interspersion 

3. Measurement of Juxtaposition 

4. Recognition of exclusion factors 

Basically, the four components Interact in the following way. Suitable 
land cover vegetation naps are either obtained from existing sources or 
compiled. The necessary/ vegetation categories, nap scales and ninimuo 
napping unit size nay vary from region to region and with the species for 
which potential habitat is being assessed. The habitat criteria far the 
species of interest must be known (or estimated). Such criteria include 
t.he relative desirability (i.e., the weighting factors) of various 
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v«g«cacion/liDJ cover edges sod che snisels' preference for various 
vegecaclon discribuclon paccoros. Clasalflcaclon of vegtcacion groups is 
sooier.laes racher arbitrary (Plelou, 1977), and susc be made biologically 
in terms of the requirements of the organisms involved. For example, the 
timber type classification system used by forest industries may not always 
be adequate for use in wildlife habitat inventories. Finally, specific 
restrictive factors or resources (e.g., water) that either MUST or 
MUST MOT be present for suitable habitat need to be known. 

A spatial diversity "SD" index value is computed for each parcel of land 
(cell) (of any desired size) relative to each wildlife species or groups 
of species. The index is a function of ”15," interspersion, "JX," 
juxtaposition, and any number of restrictive factors. 


^°A- 

M ^ Q t| tf 

^ A 3 ^ ^^A i: ^ 


LlI * 

0 X 

Q 



1 

1 

t 

1 ^ 

A 

A 


where: 

A « indicates a specific wildlife species or group of species (3, 

C, D, etc. for others). 

7 * indicates the relative importance of interspersion to juxtaposition 
for wildlife species A, 3, or C, etc. 

3 ■* indicates the relative importance of juxtaposition to interspersion 
for wildlife species A, 3, or C, etc. 
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Moca Chat 3 and a can rangt bacvcan 0 and <*>1 buc atusc sua co 1.00. A low 
valua would Indicaca a vat7 undaalrabla or uniaporcanc characcarlacic and 
'^1 as var/ daslrabla or laporcanc characcariatlc. Scaling will hava co 
ba workad ouc and a sanaleivicy analysis parforaad. 



Indlcacas a rascricclva factor chat Is assanclal for wlldllfa spacias 


group "A". An txaaple of a rascricclva factor nlghc be the presence of 

vacar vlchln one aile. If chls is present (l.e., satisfying a necessity), 

than I 1 1 is given a value of 1 and has no iapact on the value of IS . 

A 

l-iowever, If there is no vacar (an absolute necessity) , chan 1 1 j is assigned 

A 

a valua of "0" and aucooacically aakas IS^ *10. In soma cases the 

restrictive factors aay be sac at inceraediace values indicating undesirable 
conditions buc not total exclusion. Values for 15^ referring to "high," 
"aadiua, " and "low" hava to ba decertilnad (categorized). 


^lecessar*/ Land Cover 3a ca 


It is assumed chat a suitable vegetation nap is available which includes the 
necessary categories of overscory and/or understor/ communities indicated. 
This must be determined for each wildlife species for which habitat is to 

be assessed. 


A. Small Area, Manual Analysis . A grid drawn on clear plastic materia 
is placed directly over the vegetation nap. The predominant vegetation 
tacegorL^ in each cell of the grid is determined, and coded directly on the 
clear plastic using a grease pencil. A key will be needed to relate the 
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letters or symbols used to the vegetation categories. 

3. Large Area, Computerized Analysis . A vegetation aap in polygon 
form is digitized (or manually coded) at any desired cell size. Individual 
cells are categorized and a file created to store the restilting data. 

Measurement of Intersoerslon 

A. Small Area, Manual Analysis . The vegetation categoir/ predominant 
in each individual cell on the clear plastic grid is compared to each of 
the immediately adjacent cells. The number of adjacent cells of another 
vegetation type are counted and that number recorded on the plastic sheet 
in the lover right hand comer. 

Consider the following rvo examples: 


Exastole I 


Examole II 


3 

: 

A j A 

3 

-si 

3 1 

A i A 


IS - 3 


3 

c 

1 

3 i 

: A 

A, 

/ 

3 ! 

i 3 

C 

I 


IS - 7 


The center cell in example 1 has 3 adjacent cells vlth dissimilar predominant 
vegetation types. Therefore, the value for interspersion is 3. In the 
second example, the IS value is 7. It is clear that the land cover 
patterns are much more intermixed in example II. Tnose ceils with IS 
values of 7 or 3 could be printed light gray, values of 3-6, intemeciata 
gray, and 0-2 as dark gray. Mote that each cell in the entire matrix 
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becomes the cencrold cell for comparison vrlch adjacent cells. Thus, a map 
of icterspersion is produced from all of the "IS" values computed by oioving 
the 3x3 matrix throughout the data set. 

3. Large .\rea. Computerized Analysis , k computer could easily be , 
prograsnned to compare adjaceut cells and create a file vlth the interspersion 
values for each cell. Any range of IS values cotald be assigned a specific 
color or gray tone, and thus an Interspersion map could be made. 

Measure or Juxtaposition 

Wildlife biologists know that certain types of vegetation edges are very 
important for specific wildlife species. Abundance of these species 
may be considered a consequence of edges where types of food and cover come 
together (Leopold, 1936). According to Odum (1971), the edge effect may 
be defined as the tendency for an increase in variety and density of 
organisms at community junctions. This effect is most nark^-* in animals 
with relatively law mobility (Leopold, 1936) and high requirements in 
terms of diversity of vegetative communities (Leopold, 1936). Various 
edge combinations can be assigned a relative weighting factor for each group 
of wildlife, e.g., 

A/3 .60 

A/C .30 

3/C .10 

In this case the relative value of an A/ 3 edge is twice that of an A/C 
edge for a particiular wildlife species. Therefore, a measure of juxta- 
position can be easily computed by sumaiag the various quantity-quality 
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produces for all edges relative to each centroid cell In the data aatrlx. 
Considering exaople aatrlces I and II again: 


Example I 


B 

A - A 

t 

3 1 4 

3 1 A 

A 

3 ! 

A 


Edge 

Type 

Quantity 

Quality 

Total 

A/ 3 

4* 

.60 

2.40 

A/C 

0 

.30 

0.00 

3/C 

0 

.10 

0.00 

JX Index - 


2.40 


Example II 


Edge 

Quantity 

Quality 

Total 

.V3 

5 

.60 

3.00 

A/C 

5* 

.30 

1.50 

3/C 

0 

.10 

0.00 

JX Index ■ 


4.50 


3 

1 

C 

3 

A 

A 

3 

1 

3 

I c 

c 

1 


The JX value for example I is 2.10 and 4.30 in example II which has more 
edges which are of importance to the wildlife species -under consideration. 
* Note that diagonal edges only count 1 while either vertical or horicenta 


edges count as 2. 
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RESULTS AND DISCUSSION 

A pcrcion of che original coded vegecacion aap and che resulting maps for 
incerspersion, juxtaposiclon and spatial diversity are shown in Figures , 

1, 2, 3, and 4. The area shown includes 20 rows and 28 coluans of the coded 
input data. The numbers in Figure 1 correspond to coefficients which were 
arbitrarily assigned to the various vegetation categories napped by 
Ganaon and Carter (1979). 

The dark, intemediate gray and light areas in Figure 2 represent low, 
nediua and high interspersion, respecti'/ely. These correspond to che 
following ranges for che "IS" calculation, respectively: 

0 CO .3 
>.3 CO .6 
>.6 to 1.0 

The designations of dark, intemediate gray, and light in Figure 3 show 
juxtaposition and co'^respond to these ranges for che "JX" calculation, 
respectively: 


0 CO .3 
>. 3 to .6 
>.6 CO 1.0 

Finally, che spatial diversity index "SD" was categorized in an identical 
way. The resulting nap is shown in Figure 4. 
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SUM^Y AND CONCLUSIONS 

A wildlife hablcac dlverslcy aap was produced for a hypochecical wildlife 
species in the Olsnal Swamp ucllizlag a vegetation cover aoap. This aethod 
can be performed very quickly by computer over large areas, given the 
aeccssary input data. Maps of intersperslon and Juxtaposition can be 
produced as well by assigning printer symbols to arbitrarily designated 
categories for each of the three parameters (intersperslon. Juxtaposition 
and the wildlife habitat diversity index). Such aaps are repeatable and 
would be consistent over large areas. The aost crucial part of the 
operation is the assignment of the weighting factors from "known" ecological 
information about each wildlife species. The computerized aethodolog^/ 
aay have tremendous potential when implemented with remotely sensed digital 
data far land/cover vegetation. 

Further work is needed to determine the sensitivity or the output maps to 
changes in the weighting factors for various species of wildlife. The 
relation between animal home range and suitable cell size must also be 
examined. More efficient methods should be used to digitize the land 
cover/’vegetation maps. Finally, the maps must be more thoroughly evaluated 
by field resource managers and wildlife habitat specialists. 

The method proposed here measures only the spatial diversity of the landscape. 
Such a measure, and the maps which result, could be incorporated into a 
larger, more comprehensive system for assessing wildlife habitat O'cality. 
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'isura 1. Digically ceded cellol.ar =ap o>* a sna” 
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cover aao. 
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jure 1. ".K" 

low, CIS" is >.3 zo .5), ”3" is iac3r=ediace ("SD" is 0 
3), and ”C" is high ("SD" is >.6 to 1.0) spatial divarsiry. 





















